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Blue stragglers in globular clusters are abnormally massive stars that should 
have evolved off the stellar main sequence long ago. There are two known pro- 
cesses that can create these objects: direct stellar collisions^ or binary evolution^. 
However, the relative importances of these processes have remained unclear. In 
particular, the total number of blue stragglers found in a given cluster does not 
seem to correlate with the predicted collision rate^'^, providing indirect support 
for the binary scenario. Yet the radial distributions of blue stragglers in many 
clusters are bimodal, with a dominant central peak^~^. This has been interpreted 
as an indication that collisions do dominate blue straggler production, at least 
in the high density cluster cores^'^. Here, we show that there is a clear, but 
sub-linear, correlation between the number of blue stragglers found in a cluster 
core and the total stellar mass contained within it. This is the strongest and 
most direct evidence to date that most blue stragglers, even those found in clus- 
ter cores, are the progeny of binary systems. However, we also point out that 
the parent binaries may themselves have been affected by dynamical encounters. 
This may be the key to reconciling all of the seemingly conflicting results found 
to date. 

In both blue straggler (BS) formation scenarios, the numbers and properties of the 
resultant BSs can shed light on the dynamical history of their host clusters. The conclusions 
drawn from the observed BS populations depend on which formation mechanism is dominant. 
Unfortunately, even though a number of potential formation tracers have been proposed, a 
simple observational distinction between collisional and binary BSs has yet to be found 
In general, it is impossible to ascertain whether an individual BS was formed via a stellar 
collision or a binary merger. 



However, it should be possible to determine the dominant formation channel statisti- 
cally, by studying the relationship between the size of a BS population and the properties of 
its host cluster. Here, we focus on BSs found in the cores of clusters, where collisions should 
be most frequent. If most BSs are formed as a result of single-single collisions, their number 
in a particular cluster core should be approximately Nbss,coII — tbss/tcoU, where tbss is 
the typical BS lifetime, and Tcou is the time-scale between single-single stellar collisions in 
the core. Conversely, if most core BSs are the progeny of binary stars, their number may 
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be expected to scale simply as NsssMn fbinMcore, where fhin is the binary fraction in the 
core and Mcore is the total stellar mass contained in the core. These predictions are for 
the simplest versions of the competing scenarios, but note that more complex, intermedi- 
ate channels could also contribute (e.g. BSs produced as a result of dynamical encounters 
involving binary stars). We will return to this subtlety below. 

In order to test these simple predictions, we used our existing data base of core BS 
counts which are derived from a large set of HST-hased colour-magnitude diagrams^^. 
Figure 1 confirms that there is no global correlation between the observed core BS numbers 
and those predicted by the collision scenario^' ^. In particular, the predicted numbers are 
very different for high- and low-density clusters, which is in conflict with the observations. 
By contrast. Figure 2 shows that the observed numbers do correlate strongly with the cluster 
core masses, as expected if most BSs are the progeny of binary stars. A power law fit suggests 
Nbss oc mO-=^«±o-°1 

When consideration is restricted to dense cores (po ~ 10^ Mq pc~^), the coUisional 
predictions are, in fact, significantly correlated with observed BS numbers (see the black 
points in Figure 1). However, even for these dense cores, the correlation with core mass 
remains stronger (Figure 2). Thus not only is core mass the only successful global predictor 
for core BS numbers, but it remains a better predictor than collision rate even for dense 
clusters. 

One way to distinguish even more clearly between coUisional and binary scenarios is 
to isolate the key cluster parameters that drive the observed correlations and quantify their 
impact on the observed sizes of the core BS populations. For example, even though collision 
rate and core mass are correlated parameters, only collision rate depends fundamentally on a 
cluster's central velocity dispersion. We therefore carried out a direct fit to the observed core 
BS numbers, using the generalized model Nbss — K r'^ ctq , where K, a, (3 and 7 are all 
treated as free parameters. For purely coUisional BS formation, we expect a ~ 3, /9 ~ 2 and 
7 ~ — 1, whereas pure binary formation implies a ~ 3, /5 ~ 1 and 7^0. These relationships 
assume constant binary fraction across clusters; we will return to this assumption below. 
Allowing for a more general power-law dependence on collision rate or core mass, we predict 
a ~ 1.5/3 ~ —87 for coUisional BS formation, and a ~ 3/3 with 7 ~ for the binary scenario. 

Figure 3 shows the results of this fit. The best-fit power law indices arc a = 1.47±0.29, 
f3 = 0.55 ± 0.15 and 7 = —0.40 ± 0.27. These coefficients are consistent with a power-law 
dependence on Mcore, but not with one on collision rate. (Note that since Tc, Po and a are 
correlated parameters, a, P and 7 are also correlated.) There is no evidence for a significant 
dependence of BS numbers on velocity dispersion. A fit to the dense clusters only yielded 
results which were consistent with the values obtained for the complete cluster sample. 
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We find that BS numbers depend strongly on core mass, but not on collision rate. We 
argue that this is strong evidence for a binary, as opposed to a collision, origin for blue 
stragglers in globular clusters. We note that the core mass dependence can be approximated 
as Nbss with 5 ~ 0.4 — 0.5. The sub-hnear dependence of BS numbers on core 

mass is interesting. For ordinary populations of single stars, one would expect a roughly 
linear dependence, and we have verified that this is the case for the red giant and horizontal 
branch populations in our cluster sample. In the context of the binary scenario for BSs, the 
simplest way to account for this sub-linearity is to posit that core binary fractions decrease 
with increasing core mass, i.e. fbin oc M^^^g with e ~ 0.5 — 0.6. This is a testable prediction. 
There have been two recent efforts to estimate core binary fractions in a consistent way 
from homogeneous, deep HST/ACS observations of globular clusters^^' . One of these^^ is 
restricted to a set of 13 low-density clusters (Sample 1). The other^^ is still preliminary, 
but based on a larger set of 35 clusters that span a wider range in density and other dy- 
namical parameters (Sample 2). Promisingly, from our perspective, a correlation between 
BS frequency and core binary fraction has already been reported for Sample 1^^, and an 
anti-correlation between core binary fraction and total cluster mass has been reported for 
Sample 2^^. Since Sample 1 spans only a narrow dynamic range, and the binary fractions 
reported for Sample 2 are only available in graphical form, we could only make a preliminary 
attempt to infer the dependence of fhin on Mcore from these samples. Nevertheless, our power 
law estimate is fun oc M~^f^ for both samples, close to the predicted relation, albeit with 
considerable scatter. Thus the number of BSs found in the cores of globular clusters scales 
roughly as Nbss oc fbinMcore, just as expected if most core BSs are descended from binary 
systems. 

There are several ways to build on the analysis presented above. First, binary-single 
and binary-binary interactions may actually produce more blue stragglers in clusters than 
single-single collisions because the interaction cross sections are larger. The expected num- 
ber of binary-single [binary-binary] collisions scales like the number of single-single collisions 
times fb [f^]. A better determination of the dependence of binary fractions on cluster param- 
eters is clearly warranted to test these predictions fully. However, based on the preliminary 
analysis of the observed binary fractions described above, it currrently seems unlikely that 
binary-single and binary-binary interaction numbers can produce the observed scaling of 
blue stragger numbers with cluster parameters. 

Second, the effects of mass segregation have been ignored in our analytical predictions. 
Massive stars at any particular location in a cluster will sink to the core on roughly the 
local relaxation time. Thus some of the BSs found today in cluster cores may actually have 
been formed outside the core^'®'^^. Moreover, clusters evolve, so the present-day cluster 
parameters may not always be close approximations to what should ideally be suitably- 
weighted averages over the BS lifetime. As a first test, we have split our cluster sample in 
objects with half-mass relaxation times smaller and greater than the blue straggler hfetime. 
If mass segregation were responsible for masking an underlying trend with collision rate. 
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we may expect the long relaxation timescale clusters to reveal this trend, while the short 
relaxation time clusters should exhibit enhanced BS numbers (at a given collision rate). We 
find no evidence for either of these effects. It is also worth noting that our result for core 
BSs is actually consistent with the behaviour of the total (core+peripheral) BS populations, 
which also seem to show a sub-linear correlation with total cluster mass 

Our results do not imply that dynamical interactions are unimportant to BS formation. 
The properties of the core binary populations from which BSs are probably descended may 
themselves be strongly affected by dynamical encounters. There is a wide spectrum of 
possible paths within the family of binary channels. While some 'binary BSs' may be the 
direct descendants of primordial binaries that evolved essentially in isolation, others may 
exist only as a result of multiple strong interactions and exchanges 29,30_ ^^^^ ^^le first 
to point out this subtlety (though it has sometimes been obscured by the terminology used). 
For example, Mapelli et al. note that the coUisional BS population invoked in their model 
could be produced via encounters involving binaries^' ^, and the primordial binary channel 
described by Davies et al. explicitly allows for dynamical interactions^^. Our results do, 
however, argue against the simple scaling of blue straggler numbers with any collision rate. 
Clearly, a key challenge for future work will be to determine if dynamically active binaries 
can produce BS numbers consistent with the trends found here. 
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Figure 1: The observed numbers of core blue stragglers versus the numbers expected from 
single-single collisions. Black points correspond to high-density clusters (po > 10"^ Mq pc~'^), 
grey points to low-density ones (po < 10^ pc~^). The error bars shown are approximate 
68% Poisson confidence intervals^^. There is no correlation betwen observed and predicted 
numbers across the full sample (Spearman's rank test yields ps = —0.03, p = 0.6, = 56). 
For the subset of dense clusters, a correlation does exist (ps = 0.52, p = S x 10^"^, = 34). 
In calculating the predicted numbers, we took the time-scale between collisions to be^'' 

^coii = 1.1 X 10^° {f^y {j^y" {^^) (H)"' (fe)"' years. Here r, is the core 
radius, no is the central number density in stars per cubic parsec, cxo is the central velocity 
dispersion, is the average stellar mass (taken to be 0.5 Mq), and is the average 
stellar radius (0.5 Rq). Core radii and central luminosity densities were taken from the 
McMaster Globular Cluster Catalogue^^. Core number densities were estimated by adopting 
a fixed mass-to-light ratio of 2^^. Central velocity dispersions were taken primarily from 
the compilation of Pryor & Meylan^" , but supplemented with data from Webbink^^. The 
typical BS lifetime in a cluster was taken to be 1.5 Gyr 22,23_ jsjQ^g that for high-concentration 
clusters, the observational definition of the core radius can be somewhat arbitrary. However, 
this should not have a significant impact on our analysis, as long as the associated central 
densities and velocity dispersions remain reasonable averages over the relevant regions. 
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Figure 2: The observed numbers of core blue stragglers versus the estimated core masses. 
Black points correspond to high density clusters (po > 10^ pc~^), grey points to low- 
density ones (po < 10^ Mq pc"^). The error bars shown are approximate 68% Poisson 
confidence intervals^^. BS numbers and core masses are strongly correlated (p^ = 0.71, 
j9 = 4 X 10~^°, = 56), as expected if BSs are descended from binary systems. The 
correlation with core mass also holds for the subset of dense clusters (p^ = 0.84, p = 3x 10~^°, 
= 34) and in both cases is stronger than that with collision rate. The dashed line is a 
power law fit to the full set of clusters, obtained with a method appropriate for this type of 
data set^^ The best-fitting power law is Nbss oc M^-^^^^-^^. 




Figure 3: Results of fitting the generalized model Nbss oc -fC r° pg the observed 

BS numbers. The top panels show three projections of the 4-dimensional parameter space. 
The contours in each panel represent la, 2a and 3a joint confidence intervals on the two 
parameters being considered. Black contours correspond to the full cluster sample, grey 
contours to the dense cluster sample (po > 10"^ Mq pc~^). The straight lines in the middle 
and right panels show the parameter relationships expected for a power law dependence on 
collision rate (dashed lines) and core mass (solid line). The best fit power law indices are 
a = 1.47 ±0.29, /3 = 0.55 d= 0.15 and 7 = -0.40 ±0.27. The bottom panel shows observed BS 
numbers vs the best-fit model prediction for the full sample. Black points correspond to high 
density clusters (po > 10^ pc~^), grey points to low-density ones (po < 10'^ Mq pc~^). 
The error bars shown are approximate 68% Poisson confidence intervals^^. The dashed line 
marks the identity relation, i.e. N Bss{observed) = NsssiP^^d'^'Cted) . 

We also carried out a separate fit to the subset of dense clusters (po > 10^ Mq pc~^). This 
yielded a = 1.59 ± 0.33, P = 0.63 ± 0.19 and 7 = -0.28 ± 0.31, consistent with the values 
obtained for the complete cluster sample. Note that, even within this subset, there is no 
evidence for a dependence of BS numbers on ctq. 
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